Passive scalar turbulence in high dimensions 
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Exploiting a Lagrangian strategy we present a numerical study for both perturbative and nonper- 
turbative regions of the Kraichnan advection model. The major result is the numerical assessment 
of the first-order 1/d-expansion by M. Chertkov, G. Falkovich, I. Kolokolov and V. Lebedev (Phys. 
Rev. E, 52, 4924 (1995)) for the fourth-order scalar structure function in the limit of high dimen- 
sions d's. In addition to the perturbative results, the behavior of the anomaly for the sixth-order 
structure functions vs the velocity scaling exponent, £, is investigated and the resulting behavior 
discussed. 



PACS number(s) : 47.27.Te, 47.27.-i 

A first principles theory of turbulent advection was 
very recently developed in the framework of the passive 
scalar model introduced by R.H. Kraichnan |l],|j where 
the velocity field advecting the scalar is Gaussian in space 
and white in time (i.e. ^-correlated). The main feature 
of such model is that its statistical description can be 
obtained from the solution of a set of closed linear par- 
tial differential equations for the equal-time correlation 
functions jjj. 

It was conjectured by Kraichnan g] that, despite the clo- 
sure of the equations for the equal-time correlation func- 
tions, the system should display an intermittent behav- 
ior. This means [|| that the scaling exponents £2™ of the 
2n-order scalar structure functions S2 n are anomalous, 
namely S 2n = ([6(r,t)-6(0,t)] 2n ) oc r^" with Q 2n < n( 2 
as T goes to zero. This is in sharp contrast with the nor- 
mal scaling C^n = "C2 predicted on the basis of mean- field 
dimensional arguments. 

Intermittent behavior was subsequently established by 
Gaw§dzki and Kupiaincn [0], Chertkov et al [0 and 
Shraiman and Siggia HQ- Anomalous scaling of the 
structure functions appears because of the presence of 
zero modes of the equations for the correlation functions 
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The second-order structure function was solved exactly in 
Ref. M and it does not show any anomaly. Corrections to 
normal scaling start to appear for the fourth order struc- 
ture functions (the third, when scalar fluctuations are 
sustained by a large-scale gradient The anomalous 

exponents have been up to now computed only perturba- 
tively around limit cases where the structure functions of 
the Kraichnan model are known and display normal scal- 
ing. The corresponding perturbative regions are : small 
£'s large d's |];|l(J and £ close to the Batchelor limit 
£ = 2 pj^,[l2| 13|, £ being the scaling exponent of the ad- 
vecting velocity field and d the dimension of the space. 
The first two expansions are regular, while for the third 
one the relevant small parameter should be \J2 — £. This 
is due to the preservation of the collinear geometry in the 
Batchelor limit, leading to an angular non- uniformity in 
the perturbation analysis |. 

The scaling behavior of high order structure functions 



was investigated by means of instanton calculus |lq ] . The 
solutions obtained up to now predict the saturation to a 
constant value of the scaling exponents £ n 's as n becomes 
large ||||. 

It is of great interest both to develop efficient numerical 
strategies to test the predictions of the available pertur- 
bation theories - and thus to give numerical confirma- 
tions of the mechanisms responsible for the emergence 
of intermittency - and to investigate the nonperturba- 
tive regions, still not completely accessible by analytical 
approaches. Exploiting a new numerical strategy such 
program was pursued in Refs. p9|-|2l| (see also the re- 
view paper ||22| ) where, in particular, the behavior of the 
anomaly o~2n = nC^ — ($n wa s studied as a function of £ 
for n — 2 both in two and three dimensions. The numer- 
ical experiments confirmed the perturbative predictions 
for £ — > and £ — > 2. However, due to the difficulty 
of the numerical simulations no results were available for 
the behavior of 02n as a function of £ for moments higher 
than the fourth and far from the two above perturbative 
regions. A numerical study of the curve a& vs £ in the 
nonperturbative region of the Kraichnan model is one of 
the aims of the present Letter. 

The second aim is the investigation of the Kraichnan 
model in large spatial dimensions. The 1/d expansion 
received solely a partial numerical confirmation for the 
third-order structure function ]23| ] . The value of d inves- 
tigated in Ref. [^3) was not large enough to reach the per- 
turbative regime for the fourth-order structure function 
and thus to verify the perturbative prediction of Ref. [|| . 
This is a second goal of the present Letter. Exploiting the 
Lagrangian strategy presented in Refs. [ ^9pC| ] we present 
two sets of numerical simulations which make possible the 
numerical assessment of the 1 /d-expansion of Ref. §. We 
found that, for £ = 0.8, the value of d at which the per- 
turbative regime studied in Ref. 0] takes place is d ~ 30. 
The latter value reduces as £ decreases, being of the order 
of 25 for £ = 0.6. 

Let us briefly recall the Kraichnan advection model 
Gl||. In this model the velocity field v = {v a ,a = 
1, . . . , d} advecting the scalar is incompressible, isotropic, 
Gaussian in space, white-noise in time; it has homogc- 
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neous increments with power-law spatial correlations and 
a scaling exponent £ in the range < £ < 2 : 

(K(r, t) - v a (0, 0)]Mr, t) - ^(0, 0)]) = 25{t)D a p{r), 

(1) 



where, 



(2) 



The specific form inside the squared brackets is dictated 
by the incompressibility condition on v. 
The nonintermittent velocity field defined by (|l|) advects 
a scalar field 



v • 80 = nd 2 9 + / 



(3) 



of which we want to investigate the statistical properties 
related to intermittency. Here / is an external forcing 
and K is the molecular diffusivity. 

The forcing term permits to attain a stationary state de- 
fined by the balance between production of scalar vari- 
ance (related to /) and its dissipation (related to nd 2 9). 
The net result is that (0 2 ) is finite in the steady state. 
We shall assume the (random) forcing / of zero mean, 
isotropic, Gaussian, white-noise in time and homoge- 
neous. Its correlation is specified by: 



(/(r J *)/(0,0))= X (r/£) <*(*), 



(4) 



with x(0) > and X ( r /L) a rapidly decreasing function 
for r ^> L where L is the (forcing) integral scale. 

Rather than integrating the partial stochastic differen- 
tial equation (|3|) (that is a prohibitive task already in two 
dimensions due to the (5-correlation of v) , we exploit the 
Lagrangian formulation of the passive scalar dynamics as 
m Refs. gf,§0| in order to simulate particle trajectories 
by Monte-Carlo methods. 

We recall that when such method is adopted the evalua- 
tion of structure functions is reduced to the study of the 
statistical properties of the random variable describing 
the time spent by pairs of particle with mutual distance 
less than the integral scale L. Generally, the distance be- 
tween pairs of particles tends to increase with the elapsed 
time but, occasionally, particles may come very close and 
stay so; the phenomenon is the source of scaling anoma- 
lies. The main advantage of the Lagrangian strategy is 
that one does not have to generate the whole velocity 
field, but just the velocity field over the Lagrangian tra- 
jectories. Simulations in very high space dimensions as 
well as the investigation of high-order moments become 
thus feasible. 

We now present results for structure functions up to 
sixth order. The three-dimensional results relative to the 
fourth-order structure functions have been already pub- 
lished in Ref . |nj . They will be reported in the following 
for comparison with our new data. 

The L-dependence of Se(r;L) is shown for the three- 
dimensional case in Fig. for £ = 0.9 (above) and £ = 1.5 



(below). Similar plots have been obtained also for other 
values of £. In both cases the scaling region is indicated 
by a dashed straight line the slope of which yields the 
anomaly og = 3C2 — Ce- The measured slopes for £ = 0.9 
and£ = 1.5 are cr 6 = 0.91 ±0.03 and cr 6 = 0.58 ±0.02, re- 
spectively. The error bars are obtained by analyzing the 
fluctuations of local scaling exponents over smaller ratios 
of values for L. As explained in Ref. the number of 
realizations needed to obtain high quality scaling as the 
ones displayed in (l9) increases rapidly with decreasing 
£. This is the reason why we restricted our analysis to 
the interval £ > 0.5. 

The results for the anomalies vs £ are summarized in 
Fig. §. The upper curve refers to a^, lower curve to 04 
(the latter has been already published in Ref. |]l9|). 
Some remarks on the behavior of ctq are worth. First, we 
note that when £ decreases from 2 to the anomaly a§ 
grows at first, achieves a maximum and finally decreases 
as in the case of 174. 
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FIG. 1. Three-dimensional sixth-order structure functions 
Sa vs L for £ = 0.9 (above) and £ = 1.5 (below). The dashed 
straight lines are the best fit slopes calculated in the scaling 
region. Their values give the anomalies reported in Fig. 
(upper curve). Number of realizations ~ 10 7 . 
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FIG. 2. The anomaly erg = 3(^2 — Cs f° r the sixth-order 
structure functions in three dimensions (upper curve) and, for 
comparison, the a 4 curve already published in Ref. JlSJ] (lower 
curve). The dashed straight line on the left is the first-order 
prediction 05 = 12£/5 by Bernard et al Notice the shift 
on the left of the maximum of the o§ curve with respect to 
the maximum of 04. 

We have evidence that ae tends to vanish for £ — > as 
follows from the perturbative predictions |||| the leading 
order of which o-q = 12£/5 is shown as a dashed straight 
line on the left of Fig. |[ The fact that the anomaly for Sq 
is higher than the one for S4 is an immediate consequence 
of Holder inequalities f| . 

More interesting is the fact that the maximum of the 
anomaly occurs for a value of £ which for Sq is smaller 
than for S4. This can be explained as follows. Near £ = 
the dynamics is dominated by the nearly ultraviolet- 
divergent eddy diffusion. The latter arises from small- 
scale fluctuations averaged over time scales much larger 
than the typical ones. Large fluctuations are thus aver- 
aged out and appear strongly depleted. On the contrary, 
near £ = 2 the dynamics is dominated by the nearly 
infrared-divergent stretching. No averaging of the large 
fluctuations now takes place and the net result is that 
they strongly contribute to the scalar statistics. Since 
higher moments are more sensitive to large fluctuations 
than the lower ones, one expects stretching to dominate 
the dynamics for wider ranges of £ as the order n of the 
moment increases. As a consequence the value of £ where 
the effects of stretching and diffusion balance (i.e. where 
the maximum of the anomaly is expected) should move 
toward the left as n increases in agreement with the re- 
sults reported in Fig. |[ 

Let us now turn to the behavior of fourth-order struc- 
ture functions for large space dimensions. In such limit 
the first-order perturbative prediction in X/d gives for the 
fourth-order structure functions the anomaly 4£/d P,^0[. 
The Lagrangian strategy allows to verify numerically the 
perturbative prediction. 
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FIG. 3. Fourth-order structure functions for £ = 0.8 and 
d ranging from 5 to 30. Dashed straight lines represent the 
anomaly from the 1 / d-expansion by Chertkov et al j^] . Num- 
ber of realizations ~ 10 6 . 
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We performed two sets of simulations for £ = 0.6 and 
£ = 0.8 and different spatial dimensions, d, from 5 to 
30. The fourth-order structure functions for d = 5, 10, 
20 and 30 are shown in Fig. || for £ = 0.8. Similar scal- 
ing laws have been observed also for £ = 0.6. Dashed 
straight lines represent the slopes 4£/<l For d = 30 the 
anomaly obtained from numerical simulations (i.e. the 
slopes of the curve with circles) is practically indistin- 
guishable from the perturbative expression. The discrep- 
ancy increases rapidly as d reduces. The results are sum- 
marized in Fig. Q where the anomaly 04 is shown as a 
function of d for £ = 0.6 and £ = 0.8. The fact that for 
£ = 0.6 the perturbative regime starts for a value of d 
(~ 25) lower than the one at £ = 0.8 (d ~ 30) is the 
consequence of the fact that the small parameter in the 
perturbation theory is oc l/[d(2 — £)] rather than 1/d, 



sixth-order structure functions vs the velocity scaling ex- 
ponent £. We identified and discussed two competing 
mechanisms which control the position of the maximum 
of the anomaly along the £-axis. 
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tute for Theoretical Physics of the University of Califor- 
nia, S. Barbara. Simulations were performed within the 
INFM Parallel Computing Initiative. 
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FIG. 4. The anomalies 04 vs d for £ = 0.6 and £ = 0.8. 
Dashed lines are the anomalies 4£/d from the O(lfd) pertur- 
bation theory of Ref. [^| . 

In conclusion, we verified the 1/cZ-expansion for the 
fourth-order structure functions. We presented two set 
of simulations corresponding to two different values of 
the velocity scaling exponent £. For £ = 0.8 the pertur- 
bative regime sets in for d ~ 30. This value reduces at 
d ~ 25 for £ = 0.6. The result is expected as l/[d(2 - £)] 
is the relevant small parameter in the perturbation the- 
ory. We also studied the behavior of the anomaly for the 
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